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ABSTRACT 



The influence of broken structural homology, and the implied broken dynamical 
I ^ ■ homology, is examined for the Fundamental Plane (FP). Requiring a symmetrical 

treatment of the FP variables, a bisector method of linear regression was applied, 
in 3-dimensions, to derive the best FP. A bootstrapping procedure has been used 
to estimate the uncertainties associated with the slope of the FP. For 25 E and SO 
! j Virgo galaxies, the 'standard' FP, constructed using R X I A model parameters for the 

effective radii (-R e ,4) and the mean surface brightness within this radius (S e ,4) and 
^•f-} ' using central velocity dispersion (CVD) measurements (oo), gave a relation described 

! by #0,4 OC o-l-10i0-14g-0.56i0.09_ Usmg Sergic R l/n Ught profile model paramet ers 

and the projected, infinite aperture, velocity dispersion (atot.n), derived from applica- 
' _ ' , tion of the Jeans equation to the observed intensity profiles, gave an 'improved' FP 

described by the relation i? c , n oc cr t 1 o ? 7 n ±016 E-°- 76±0 05 . This result, based on indepen- 
| dent data, supports the previous finding by Graham & Colless (1997a) that assump- 

tions of structural and dynamical homology are partly responsible for the departure 
*-£h ' of the observed FP from the plane expected by the virial theorem, which predicts 

R oc cr 2 E _1 . Upon removal of the known SO galaxies from the sample of Virgo galax- 
Q ■ ies, the above planes were observed to change to i? c ,4 oc ctq' 19±0 ' 21 £~4' 60 11 and 

+3 , i? e ,n CT t 1 o [ 2 n ±0 ' 24 S~° 74±0 09 . The perpendicular rms residuals about these planes are 

0.084 and 0.050 dex, respectively, 
j The Fornax cluster was similarly treated, although removal of the SO galaxies left 

£> . a sample of only 7 ellipticals which had published CVD measurements. Treating the 

' range of structural and dynamical profiles present in this sample produced a FP given 

by the relation i? e>n oc a to °„ °' 78 E~ 1 1 - CI7±0 - 30 , in tantalizing agreement with the plane 
■ expected from the virial theorem, but with discouragingly large errors due to the small 

sample size. Similarly to Virgo, the perpendicular rms residuals about this plane is 
0.050 dex. 

The FP was also constructed with the purpose of using it as a distance indicator, 
achieved by minimising the distance-dependent quantity logi? against the distance- 
independent quantities £ and log a. A Virgo- Fornax distance modulus was computed 
using Working- Hotelling confidence bands (Feigelson & Babu 1992). The 'standard' 
FP parameters gave a value of 0.45±0.16 mag, where-as the 'improved' FP parameters 
gave a value of 0.25±0.12 mag. However, a full treatment of the uncertainties on the 
FP slopes, derived through a bootstrapping procedure of the 3-dimensional FP data 
set, revealed that the analytical expressions for the uncertainties on the estimated 
distance moduli, given above, should be increased by a factor of ~5. 

Key words: galaxies: fundamental parameters - galaxies: structure - galaxies: kine- 
matics and dynamics, galaxies: elliptical and lenticular, cD - distance scale 



1 INTRODUCTION 

Our knowledge about the structure of elliptical galaxies has 
progressed considerably in the last decade. No longer are 
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we satisfied with fitting a de Vaucouleurs (1948,1953) R 1 ^ 4 
intensity distribution to the light profiles of elliptical galax- 
ies. Modern imaging techniques are providing ever-more pre- 
cise measurements, revealing a range of light profile shapes 
which have been shown to be better described by the Ser- 
sic (1968) R 1/n model (Capaccioli 1989; Caon, Capaccioli & 
D'Onofrio 1993; Graham et al. 1996, and references within). 
Through varying the parameter n in the Sersic model, one 
can reproduce an exponential light distribution (n=l), the 
classic de Vaucouleurs profile (n=4), approximate a power- 
law for large values of n (approximately >15), or any of the 
intermediate forms which galaxies exhibit (Michard 1985; 
Schombert 1986; Binggeli & Cameron 1991; Caon et al. 1993; 
Graham & Colless 1997a, hereafter Paper I, to mention a 
few). Consequently, for those galaxies which are not well 
described by an R 1 ^ 4 profile, the effective radius and mean 
surface brightness within this radius will be different when 
obtained using the de Vaucouleurs model and the Sersic 
model. More importantly, these differences have been shown 
to have a systematic trend: galaxies with n<4 have their ef- 
fective half-light radii over-estimated by the R 1 / 4 model, 
and galaxies better described with values of n>4 have their 
effective half-light radii under-estimated by the R 1 ^ 4 model. 
Similarly, the mean surface brightness is under- and over- 
estimated by the 7? 1//4 model when n<4 and n>4 respec- 
tively (Graham & Colless 1997b). 

Observed variations in the galaxy structure, as de- 
scribed by the R 1 ^" model, imply a variety of dynami- 
cal structures (Ciotti 1991; Paper I) are required to main- 
tain these predominantly pressure-supported stellar sys- 
tems. This conjecture is borne out in the observations of 
elliptical galaxy dynamics (Davies et al. 1983; Illingworth 
1983; Capaccioli & Longo 1994) and has implications for the 
construction of the Fundamental Plane (FP). In particular, 
the use of a galaxy's central velocity dispersion (CVD, cro) 
to represent the overall kinetic energy of a galaxy has been 
shown to be inappropriate, (J0rgensen, Franx & Kjasrgaard 
1993; Bender, Saglia & Gerhard 1994). The velocity disper- 
sion derived from the kinetic energy of random motion over 
the entire galaxy surface has been shown to be proportional 
to erg' 6 rather than a% (Busarello et al. 1997). A procedure 
is described in Paper I, based upon the theoretical study of 
Ciotti (1991), which enables the velocity dispersion to be 
measured in a consistent fashion between galaxies. Rather 
than simply using CVD measurements, that sample differ- 
ent fractions of different galaxies' velocity dispersion profiles 
(in terms of the half-light radii) , through application of the 
Jeans equation, Paper I presented a method to compute the 
asymptotic value of the velocity dispersion within a pro- 
jected aperture of infinite radius. This measure, denoted by 
a tot, has the added significance of being equal to one-third 
of the virial velocity for spherical systems (Ciotti 1994). 

In Paper I the light profile data (aperture magnitudes 
from Bower, Lucey & Ellis 1992) from a sample of 26 
E/SO Virgo galaxies was analysed. Departures from the R 1 ^ 4 
model were found to be well-described by the Sersic model, 
and the implied non-homology in the velocity dispersion pro- 
file was explored. Using the Sersic model parameters, rather 
than R 1 / 4 model parameters, and using the computed, infi- 
nite aperture velocity dispersion rather than the CVD, the 
FP was observed to change from fl^cx^- 33 * ' 1 '^"- 79 * ' 11 
to R c , n ocal t 4 ^ '"S-f 3 "'" 8 . Broken structuraland dy- 



namical homology appeared to be responsible for some of 
the tilt of the FP, i.e. its departure from the plane predicted 
by the virial theorem, i?occr 2 E _1 . 

Ultimately, this broken dynamical and structural ho- 
mology needs to be dealt with if one is to understand the 
nature of elliptical galaxies and use them for measuring dis- 
tances (Dressier et al. 1987). Given the substantial improve- 
ments in describing elliptical galaxies, it is apt that the FP, 
and its application to distance estimates, is revisited. Us- 
ing a new data base (taken from Caon et al. 1994) obtained 
by different observers with a different instrument and us- 
ing a different reduction procedure, the influence of broken 
structural and broken dynamical homology upon the FP is 
again investigated. The data used in this paper is the surface 
brightness profiles along the major axis of each galaxy, as 
opposed to the integrated aperture magnitudes used in Pa- 
per I. The influence of SO galaxies upon the slope of the FP 
is also investigated, with comparison of the FPs constructed 
with, and without, their inclusion. 

This paper presents the FP for the galaxy clusters For- 
nax and Virgo, two key clusters in the extra-galactic dis- 
tance scale (Jacoby et al. 1992; McMillan, Ciardullo & Ja- 
coby 1993; Bureau, Mould & Staveley-Smith 1996). These 
clusters are sufficiently distant to possess significant mo- 
tion due to the Hubble expansion and yet they are close 
enough that their distances may be obtained with a plethora 
of independent observational techniques. The following sec- 
tion presents the galaxy sample and the observational/model 
data which is used in this study. The Fundamental Plane is 
computed, and its method of construction, including a new 
error analysis routine, is described in Section 3. Using two 
different methods of analysis, a distance modulus between 
Virgo and Fornax is calculated in Section 4. This result is 
compared with other estimates of the Virgo-Fornax distance 
modulus in Section 5, along with a discussion of the other 
results. Section 6 provides a summary and concluding re- 
marks. 



2 GALAXY DATA 
2.1 The galaxy sample 

The galaxy sample from Caon et al. (1994) has been used for 
this investigation. This data set has been used in a bench- 
mark study (Caon et al. 1993) which gave insight into the 
nature of elliptical galaxy light profiles by quantifying their 
deviations from the R}l 4 model. Knowing that this galaxy 
sample does not posses structural homology amongst its 
members makes it an ideal sample to explore the necessar- 
ily implied broken dynamical homology and the combined 
influence of these effects upon the FP. From the sub-sample 
of galaxies for which major-axis profiles have been modeled 
(see Caon et al. 1993 for a discussion of some of the prob- 
lems in doing this), those which have no published velocity 
dispersion data have been excluded here. This reduced the 
sample of Virgo galaxies from 33 (Caon et al. 1993) to 25, 
and the sample of Fornax galaxies from 16 (D'Onofrio, Ca- 
paccioli & Caon 1994) to 10. The 10 Fornax galaxies were 
further reduced to a sample of 9 after the exclusion of NGC 
1399, which, with its shallow light profile (n>15), could not 
have its velocity dispersion profile accurately modeled due 
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to the large uncertainties associated with the effective radius 
R e . The remaining galaxies are listed in Table 1. 



2.2 Surface brightness profiles and structural 
model parameters 

The B band surface brightness profiles from Caon et al. 
(1994) have been used in this study. They were constructed 
by these authors using the 'global mapping' procedure of 
Capaccioli & Caon (1989) which couples CCD data for the 
bright inner galaxy profile with photographic data for the 
outer profile. This technique allows very deep photometry 
of galaxies, and alleviates the problem of uncertain sky sub- 
traction errors that is associated with the use of small area 
CCD's. The method gives internal errors of less than 0.1 
mag at /j,b=26 mag. 

In this study the major-axis surface brightness profiles 
have been truncated at the same inner and outer limits that 
were applied by Caon et al. (1993) and D'Onofrio et al. 
(1994) for the Virgo and Fornax galaxies respectively. In 
general, this meant fitting the entire light profile after the 
exclusion of the inner portion affected by seeing, and the 
outer portion where sky subtraction uncertainties became 
large (>0.25 mag). This left atypical range in surface bright- 
ness from jj,B =20 to 27 mag. 

The Sersic (1968) R 1/n light profile model was fitted to 
the remaining profile, with the intensity / given as a function 
of radius R such that 



I(R) = I e exp 



-b(n) 



l/r, 



- 1 



(1) 



I e is the intensity at the radius R e , and b(n) (~1.9992n- 
0.3271, Capaccioli 1989) is defined such that R e is the ra- 
dius that effectively encloses half of the total light of the 
profile model (Caon et al. 1993; Graham & Colless 1996). 
The 'shape parameter' n describes the curvature of each 
light profile. Application of a non- linear least-squares fit of 
the R}J n model to the galaxy surface brightness data, with 
equal weight for each data point, produced similar values for 
the shape parameters to those derived by Caon et al. (1993) 
and D'Onofrio et al. (1994), confirming that the models are 
being fitted consistently between authors. 

Table 1 lists the parameters from the best fitting R 1 ^ 1 
model and also those from the best fitting _R 1//4 model, which 
was fitted in an identical fashion. E e is the mean surface 
brightness within the effective radius 7? e , computed using 
the technique described in Appendix A of Paper I. The sub- 
scripts n and 4 on the model parameters are used to denote 
the model from which the particular parameters have come; 
n signifying R 1 ^" model parameters, and 4 signifying R 1 ^ 4 
model parameters. 

Ciotti (1991), and also Paper I, explored the relation- 
ship between the projected and spatial (deprojected) quan- 
tities. By definition, the surface brightness terms had a one- 
to-one mapping, while there was an approximately constant 
offset (log.R~0.13, Ciotti 1991) between the effective radii. 
Therefore either of these quantities could be used for the 
construction of the FP. The deprojected velocity dispersion 
at the effective radius was found not to have such a simple 
relation with its projected counterpart, the central veloc- 
ity dispersion measurement. This deprojected velocity dis- 
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Figure 1. The shape parameter n is plotted against the R 1 /" 
model effective radius R e for the Virgo and Fornax galaxies listed 
in Table [l]. The associated 99% confidence interval is also drawn, 
as detailed in the text. 



persion is preferred over the CVD measurement because it 
not only allows for the range of velocity dispersion profiles 
present, but also samples different galaxies in a consistent 
fashion (CVD measurements sample typically anywhere be- 
tween 0.01-0.2 R e ). However, due to the physical significance 
of the projected velocity dispersion within an aperture of in- 
finite radius, this term, which also samples different galaxies 
in a consistent fashion, is preferred in the analysis which fol- 
lows. 

Following Caon et al. (1993) (their Fig. 3), the light pro- 
file shape parameter is plotted against galaxy scale-size in 
Figure l]. Although our methods of parameterisation differ 
(Caon et al. 1993 have used model-independent estimates 
for the half-light radius), the same trend of increasing val- 
ues of n with 7i e is evident here. The present study also 
includes Fornax galaxies as well as those from Virgo. Also 
shown in the figure are the 99% confidence intervals for the 
model parameters. These were computed by evaluating the 
chi-square values for a fine grid of points around the opti- 
mum solution. The contours drawn correspond to Ax 2 =9.21 
(i.e. Act), scaled after normalisation of the reduced x 2 of the 
optimum solution. Q 

In exploring the sensitivity of the it 1 '™ model parame- 
ters to the radial range of the surface brightness profile, the 
Sersic model was fitted to that data within the inner 40". 
This outer truncation is the same as the outermost data 
points used in Paper I, enabling a further comparison of the 
present data with that used in Paper I. The inner trunca- 
tion used by Caon et al. (1993) and D'Onofrio et al. (1994) 
was again employed, except for NGC 4473 and NGC 4636 
for which only the inner 3" was excluded. In Figure ^, the 
resulting model parameters are shown plotted against those 
obtained using the entire surface brightness profile (Table 



t It is noted that these confidence intervals only encompass the 
uncertainty of the parameters derived from the analytical fit to 
the data. They do not incorporate measurement errors or zero 
point offsets which may be present in the data. 



© 0000 RAS, MNRAS 000, 000-000 



4 A . W. Graham 



n 



(a) 
5 



10 



10 



o 1 
ft 



£0.5 
K 

£ oh 



-0.5 



—i > i > > i 

NGC 4552 ▲ 



■ • 



■ • 



NGC 4339 



▲ . NGC 4550 K 



★ 

1 1 1 'l 



NGC 4550 ^ NGC 13B0 



n — i — n — r — i i i | ■ 



NGC 4361 s ▲ 



NGC 1380. 



NGC 4339 



-0.5 



0.5 

log R / \ 
to e(a) 



1.5 



(kpc) 



Figure 2. The Sersic il 1 ' ™ model parameters from the fit to the 
entire light profile (subscript b) are plotted against those param- 
eters obtained from fit to the light profile within 40" (subscript 
a). The circles represent Virgo E's (18), the triangles are Virgo 
SO's (7), the squares are Fornax E's (7) and the stars represent 
Fornax SO's (2). 



1). It is noted that the light profiles from the smaller radial 
extension (<40") may not be representative of the galaxy as 
a whole, and therefore may not represent the over-all global 
galaxy shape. This is more likely to be true for the larger 
galaxies, especially those having _R e >40". It is also true that 
the Sersic index n is less sensitive at larger values (Graham 
et al. 1996). Thus, Figure ti offers an alternative estimate 
for the uncertainty of the R model parameters. Some of 
the more discrepant points have been labeled in the plot. 
It is noted that they are predominantly the SO galaxies. If 
one is to think of these galaxies as possessing an R 1 ^ 4 (n=4) 
central bulge and an exponential (n=l) outer disk, or some 
approximately equivalent Sersic form, then this could ac- 
count for the larger index obtained using the smaller radial 
range for the SO galaxies NGC 1380, 4339, and 4550. In fact, 
such measurements offer themselves as possible diagnostics 
for distinguishing SO galaxies from ellipticals. 

For the Fornax galaxies, the tabulated extinctions range 
from y4b=-0.05 to -0.09. Some of the Virgo galaxies also 



have negative extinctions, with a mean Virgo cluster value 
of As=0.04±0.07 (standard deviation) for the sample. The 
negative values are due to the variable ratio of Hj gas-to- 
dust in the Galaxy (Burstein & Heiles 1978), and have no 
physical meaning. Such values arise in the method used by 
Burstein & Heiles for regions where there is zero or negligible 
reddening, and they should be set equal to 0.00. Corrections 
for Galactic extinction (Burstein & Heiles 1984; Burstein 
et al. 1987) were not applied by Caon et al. (1994). This 
practice has been maintained here. 

At a redshift of around 0.003, the k-correction and red- 
shift dimming will be equally small for both clusters, at 0.015 
mag and 0.013 respectively. Consequently these terms have 
not been incorporated into the analysis. 

2.3 Dynamical model parameters 

CVDs are typically used as a measure of the kinetic energy 
within elliptical and SO galaxies. There are many problems 
associated with doing this. Firstly, rotational energies, which 
are likely to be significant for SO galaxies, are ignored. Sec- 
ondly, fixed aperture sizes (on the sky) in which the CVD is 
measured correspond to different physical sizes for galaxies 
in clusters at different distances. This effect alone may likely 
introduce systematic errors into cluster distance determina- 
tions which have incorporated fixed aperture sizes to mea- 
sure the CVD. Thirdly, given that galaxies have a range of 
sizes and structural profiles, they will also have a range of 
dynamical profiles which support the various structures that 
are observed. Therefore, even within clusters, fixed aperture 
sizes will sample different fractions of individual galaxies, as 
measured in terms of their effective radii. 

In an attempt to find a better measure of the kinetic 
energy within each galaxy, a method of computing the ve- 
locity dispersion in a consistent manner for all galaxies was 
presented in Paper I. Deprojecting the galaxy light profile, 
the spatial luminosity density (Ciotti 1991) was computed. 
Assuming a constant M/L ratio within each galaxy, such 
knowledge of the internal density profile coupled with the 
application of the Jeans hydrodynamical equation enables 
one to compute the internal velocity dispersion profile. This 
can then be projected outward to give the projected radial 
velocity dispersion profile, which in turn can be integrated 
over increasing aperture sizes to give the projected aperture 
velocity dispersion profile. This is then calibrated using the 
measured CVD. 

The computed spatial (deprojected) velocity dispersion 
at the spatial half-light radius was used in Paper I to repre- 
sent the kinetic energy of the galaxies. Unlike use of CVD, 
this approach consistently samples the velocity dispersion 
profile in the same way for all galaxies. As an alternative, 
the computed velocity dispersion within an aperture of infi- 
nite radius was also used. This quantity has been shown by 
Ciotti (1994) to be equal to one-third of the virial velocity 
dispersion, independent of any orbital anisotropies that may 
be present. This measure for the kinetic energy is adopted 
here, not only for the consistent fashion in which galaxies 
are sampled but also due to the significance of the infinite 
aperture velocity dispersion a to t- These values are listed in 
Table 1, where the subscripts 4 and n refer to use of the 
R 1 / 4 and R 1 ^" models respectively. Also listed are the CVD 
measurements, ao, from McElroy (1995). 



© 0000 RAS, MNRAS 000, 000-000 



A Virgo-Fornax distance modulus 5 



Table 1. Galaxy model parameters. The galaxy sample is a selection from Caon et 
al. (1994). The parameters are from the best fitting i? 1 / 4 and R 1 /™ models, with 
the exception of <jq which is the tabulated central velocity dispersion from McElroy 
(1995). R eA and R e , n are the effective radii of the R}^ and the R^^ n models respec- 
tively, with £ e ,4 and E e , n the mean surface brightness within R e . (r tot! 4 and <Ttot,n 
are the infinite aperture velocity dispersions described in the text. 
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3 CONSTRUCTION OF THE FUNDAMENTAL 
PLANE 

The quantities which are used in this investigation of the FP 
are the effective radius, R e , and the mean surface brightness 
within this radius, E e , derived from the model fits to the sur- 
face brightness profiles. These values are listed in Table 1. 
This approach differs from D'Onofrio, Longo & Capaccioli 
(1996), who used model-independent values for these quan- 
tities in their construction of the FP. They also only used 
CVDs for their velocity dispersion measure. 

In this section, the Fundamental Plane is constructed 
with the goal of comparison with theoretical predictions. 
Consequently, the 3 parameters which are used to define the 
FP (cr,7?,E) are treated as independent of each other. This 
requires a symmetrical treatment of these parameters, and 



the bisector method of linear regression is the preferred tech- 
nique (Isobe et al. 1990; Feigelson & Babu 1992). This ap- 
proach requires that three FPs be computed, each of which 
is obtained by minimising the residuals of a different param- 
eter. The mean slope angles of these 3 planes are then used 
to represent the final FP. 

The uncertainties associated with the FP slope are also 
important. Feigelson & Babu (1992) not only show that the 
standard formula (e.g. Bevington 1969) for calculating the 
uncertainties on the derived slope and intercept is in general 
mathematically incorrect, but they also show that the de- 
rived asymptotic, analytic expressions for the variance will 
under-estimate the true variance when the sample size is 
small (approximately <50). They therefore recommend the 
use of a bootstrapping procedure to gauge the variance on 
the fitted slope and its intercept. This approach was followed 
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Table 2. Fundamental Plane _Roccr A E s for a sample of 25 E &: 
SO Virgo galaxies constructed using the bisector method of linear 
regression in 3 dimensions. The errors are from a bootstrapping 
procedure. The approach referred to as K 1 '^[(To] has used the ef- 
fective radius from the R 1 / 4, model, the mean surface-brightness 
within this radius, and the projected CVD, cro, in the construc- 
tion of the FP. R 1 ' n [er totjn ] uses the infinite aperture velocity 
dispersion term, <r to t,n, and the effective radius and associated 
surface brightness term from the R}/ n model. 



model 


Rota A 




_L rms 




A 


B 


residual 




1.10±0.14 


-0.55±0.09 


0.099 




1.12±0.17 


-0.64±0.07 


0.084 




1.17±Q.15 


-0.60±0.08 


0.100 


# 1/n ,crtot,n 


1.37±0.16 


-0.76±0.05 


0.077 



Table 3. Fundamental Plane Rcca A H B for a sample of 18 E 
Virgo galaxies constructed using the bisector method of linear 
regression in 3 dimensions. The errors are from a bootstrapping 
procedure. 



model 


Roca 




_L rms 




A 


B 


residual 


RV\*o 


1.19±0.21 


-0.60±0.11 


0.084 


R 1/n ,cr 


1.22±0.25 


-0.66±0.12 


0.071 


R 1/4 ,VtotA 


1.28±0.24 


-0.65±0.11 


0.085 


R 1/n ,<Ttot,n 


1.72±0.24 


-0.74±0.09 


0.050 



Table 4. Fundamental Plane Raccr A T, B for a sample of 7 E For- 
nax galaxies constructed using the bisector method of linear re- 
gression in 3 dimensions. The errors are from a bootstrapping 
procedure. 



model 


Roca A 




_L rms 




A 


B 


residual 




1.69±0.83 


-0.85±0.30 


0.064 


R 1/n ,<ro 


1.79±0.79 


-0.98±0.28 


0.049 




1.85±0.86 


-1.00±0.35 


0.070 


fl 1/n ,<Xtot,n 


2.03±0.78 


-1.07±0.30 


0.050 



Table 5. Fractional Variance from the Principal Component 
Analysis on the sample of 25 Virgo E & SO galaxies. 



model 


Vari 


Var2 


Var3 




76.03% 


21.99% 


1.98% 


R 1/n ,cr 


82.79% 


15.80% 


1.41% 


i? 1/4 ,(Ttot,4 


72.30% 


25.57% 


2.13% 


R}/ n , Otot.n 


75.18% 


23.44% 


1.38% 



the data sets are defined by a 2-dimensional plane within 
the 3-space of observables. Table [| displays the fractional 
variance of the data along the major eigenvectors of the 
FP parameter space. This reveals that greater than 97% of 
the variance in the data indeed lies in a plane, confirming 
the appropriateness of constructing a plane to describe the 
properties of elliptical galaxies. 



in this study. However, this numerical resampling technique 
was modified for application to the small sample of Fornax 
elliptical galaxies in order to exclude 10% of the 400 boot- 
strapping simulations with a negative, rather than positive, 
slope. These simulations were excluded from the computa- 
tion of the final slope variance due to their obviously wild 
departure from the expected value. Such FP's are due to the 
small sample size of 7 points in Fornax. With the somewhat 
larger sample of Virgo galaxies (25) this problem did not 
arise in any of the 2000 simulations used. 

The total sample of Virgo galaxies contains 18 E's 
and 7 SO's. It was therefore possible to compute the 
FP using both the entire sample, and using only el- 
liptical galaxies. In addition, given the four sets of 
FP parameters (log a , log Re,4, E e , 4 ), (log a , log J? c , n , S c , n ), 
(logfj t ot,4,logi? c ,4, E c , 4 ) and (log cr tot , n , log i? c , n , E c , n ), a to- 
tal of eight FP's have been computed for the Virgo clus- 
ter. The results are shown in Table ^| and Table S where 
the rms residual is that perpendicular to the fitted plane. 
The uncertainties on the estimated slopes are those from 
the bootstrapping procedure, and are therefore larger than 
the standard analytically determined uncertainties derived 
for the FP. Using the sample of 7 available Fornax elliptical 
galaxies an additional four planes were computed using the 
above mentioned four sets of parameters. Table ^ shows the 
results. 

A Principal Component Analysis (PCA) was performed 
on each of the four 3 dimensional data sets from the com- 
bined E fc SO Virgo galaxy sample. The code from Murtagh 
& Heck ( 1987 ) was implemented to show the degree to which 



4 THE VIRGO-FORNAX 
DISTANCE-MODULUS 

Only the elliptical galaxy population is used in the following 
determination of the Virgo-Fornax distance-modulus due to 
the contaminating influence of the SO galaxies on the FP. 
Contaminating in the sense that SO galaxies may possess a 
significant degree of rotational energy which is ignored in 
the present analysis. 

The problem of computing the offset between two par- 
allel data sets, with an emphasis on cosmic distance scale 
calibrations, is well described in Isobe et al. (1990) and 
Feigelson & Babu (1992). If the objective in construct- 
ing the FP is for the purpose of computing distances be- 
tween galaxy clusters, then the FP parameters should be 
treated differently than they were treated in Section ^ 
(where the FP was constructed for comparison with a the- 
oretical prediction). In this case, one wishes to predict the 
value of a distance-dependent quantity, namely log R, from 
the distance-independent quantities log a and E. Conse- 
quently, the plane that is fitted to the data should be com- 
puted through an ordinary least-squares minimisation of the 
distance-dependent quantity, rather than using the bisec- 
tor method of linear regression which treats all variables 
equally (Section ^). This approach will minimise the un- 
certainty on the estimate of the relative distance between 
clusters. Such an approach was taken with the two ex- 
treme Virgo and Fornax data sets, (log (Jo, log J? G ,4,E e ,4) and 
(log crtot,n,log i?e,n,Ee,n)- (These are extreme in the sense of 
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Table 6. Fundamental Plane Rona A Y, B using only elliptical 
galaxies. Constructed using ordinary linear regression, in 3- 
dimensions, that minimised the log R residuals. The errors are 
from a bootstrapping procedure. 



model 


Roccr 




(log R) rms 




A 


B 


residual 


Virgo: R^/ A ,ao 


1.11±0.27 


-0.59±0.13 


0.134 


Fornax: flV4 )Cro 


1.14±0.84 


-0.59±0.29 


0.120 


Virgo: flVn.CTtot.n 


1.62±0.27 


-0.74±0.12 


0.100 


Fornax: R l / n , cr to t,n 


1.47±0.82 


-0.81±0.34 


0.110 



Table 8. Distance moduli given by 5(Alogi?), where Alogi? 
is the computed intercept offset between the cluster FPs, as de- 
scribed in the text. These have been obtained using a general- 
isation of the Working-Hotelling confidence bands (Feigelson &: 
Babu 1992), using both Virgo and Fornax as the calibration sam- 
ple. The uncertainties are from the asymptotic formula fits given 
in Tabic |. 



Calibrator 




R 1/n Wtot,n] 


sample 


mag. offset 


mag. offset 


Virgo: 


-0.46±0.23 


-0.32±0.22 


Fornax: 


+0.45±0.16 


+0.25±0.12 



the 'standard' versus 'improved' treatment of galaxy struc- 
ture and dynamics.) The best fitting planes, under this min- 
imisation scheme, are given in Table n, as are the rms resid- 
uals of the data projected along the logi? axis. 

In order to use the methods in Feigelson & Babu 
(1992), these 3-dimensional data sets have been reduced to 
2-dimensions by combining the distance-independent vari- 
ables logo - and S. These were combined in such a way as 
to preserve the previously obtained optimal linear regres- 
sion solution. This meant computing the 'mixing' value b=- 
B/(2.5A), where A and B are the FP exponents given in 
Table |^, giving the quantity log a + b<fi>, where </i> is 
the mean surface brightness in mag unitsj^] A further linear 
regression calculation was performed, minimising the log R 
residuals at the expense of the (logo- + b<n>) residuals. 
The code SLOPES, discussed in Isobe et al. (1990), was 
used to do this, and the results are presented in Table M. As 
expected, due to the choice of the mixing value, the asymp- 
totic formula gave the same slope as derived with the pre- 
vious treatment of the FP variables in 3-dimensions. The 
slopes obtained with bootstrap and jackknife sampling are 
seen to vary slightly and have different associated uncertain- 
ties, as expected for small sample sizes. The errors from the 
bootstrapping procedure for the Fornax data set have been 
restricted as discussed previously. 

The two techniques described in Feigelson & Babu 
(1992) for estimating the offset of parallel data sets have 
both been implemented here. These methods are subject to 
the condition that the Virgo and Fornax elliptical galaxy 
population define FPs which are parallel to each other. 
Given the uncertainties on the derived FP slopes (Table Q), 



S is the mean surface brightness in linear units. 



the Fornax cluster data does indeed define a plane which is 
consistent with being parallel to the Virgo cluster data. This 
is true when J? 1//4 model parameters and CVD's are used, 
and also when 7? 1// " model parameters and infinite aperture 
velocity dispersions are used to construct the FP. 

In the following method of analysis, a linear regression 
is performed on one of the cluster data sets, and the re- 
sulting regression line from this calibration is then applied 
to the second (parallel) cluster. One technique for measur- 
ing the intercept offset between parallel data sets was pi- 
oneered by Working & Hotelling (1929) and further devel- 
oped by Feigelson & Babu (1992) to permit the x-values 
(i.e. distance-independent values) to be random variables 
rather than fixed quantities. With this method, the cali- 
brated regression line for Virgo was applied individually to 
the data points in Fornax and a weighted mean offset at- 
tained; the method is discussed in Section 3.2 of Feigelson 
& Babu (1992) and their Appendix A2. This procedure en- 
compasses the fact that the confidence intervals about the 
y- value one is trying to predict, given its x- value, will be 
larger for those data points at the ends of the distribution. 
This is because the uncertainty in the slope of the calibration 
line/plane gives greater freedom of movement to the calibra- 
tion line at the ends of the distribution than in the center. 
This method of analysis gave an intercept offset, along the 
log R axis between the FP's defined by the Virgo and Fornax 
data sets, of Alogi?=-0.092±0.045 dex and -0.064±0.043 
dex for the {R}^ a ,<jo) parameters and the (R 1 ^" ,crtot,n) pa- 
rameters, respectively. 

In a similar fashion, the Fornax cluster has been used 
as the calibration sample, where the fitted regression line 
for Fornax is applied to the Virgo cluster. This was done 
with the Working-Hotelling confidence bands. Using the 
'standard' FP parameters (<To,-Re,4,£ e ,4), an intercept off- 
set of +0.091±0.031 dex was computed. Using the 'im- 
proved' FP parameters (o~tat,n,Re,n,Ee,n), an intercept offset 
of +0.051±0.025 dex was computed. These intercept offsets 
translate into a distance modulus of +0.45±0.16 mag (using 
J? 1//4 model parameters and <jq) and +0.25±0.12 mag (us- 
ing J? 1//n model parameters and atot,n)- The former estimate 
placing Fornax some (23±8)% further away than the Virgo 
cluster, compared to the latter estimate of (12±6)%. These 
solutions are fully consistent with those obtained using the 
Virgo cluster as the calibration sample (see Table ^) . 

An alternative approach for computing the offset of par- 
allel data sets was achieved following the methods in Section 
3.1 of Feigelson & Babu (1992) and their Appendix Al. The 
FP data set from the Virgo cluster of galaxies was used to 
calibrate the regression line which was then applied to the 
FP data set of the Fornax cluster as a whole (and vice-versa) . 
This was done in order to derive the intercept offset, along 
the log J? axis, between the FP's defined by the Virgo and 
Fornax clusters. The solutions obtained were similar, within 
0.001 dex, to those calculated using the Working-Hotelling 
confidence bands, although the associated uncertainties were 
generally equal or larger. This served as a useful consistency 
check on the determination of the cluster-cluster distance 
modulus. 

With the Working-Hotelling confidence bands, the un- 
certainty on the intercept offset, and hence the distance es- 
timate, is reduced by the square-root of the sample size one 
is applying the calibrated solution to. Consequently, apply- 
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Table 7. Ordinary least-squares regression analysis for logi? c versus (logcr + 6E c ), using 
the program SLOPES (Feigelson & Babu 1992). In each case, the value of b = -B/(2.5A) 
is derived from the best fitting plane shown in Table 6. 



method 


Asymptotic Formula 


Bootstrap 


Jackknife 




Slope (A) 


Slope (A) 


Slope (A) 


Virgo, RV 4 [<r ]: b=0.21 


1.11±0.08 


1.10±0.08 


1.11±0.09 


Fornax,/? 1 / 4 ^]: b=0.21 


1.14±0.27 


1.21±0.45 


1.14±0.38 


Virgo, ^/"[(Ttot.n]: b=0.18 


1.62±0.09 


1.62±0.09 


1.62±0.10 


Fornax, R 1 /" [<7 to t,n] : b=0.22 


1.47±0.38 


1.40±0.79 


1.41±0.75 



ing the Fornax regression line to Virgo, with a sample size 
of 18, results in a more accurate solution than achieved by 
applying the Virgo regression line to Fornax, with a sample 
size of only 7. However, this effect does have to compete 
with the fact that the uncertainty on the regression line is 
larger for a smaller calibration sample. As the generalised 
Working-Hotelling confidence bands propagate the errors in 
the calibrated line, this effect cannot be avoided, and the 
size of these errors influence the associated errors for the in- 
tercept offset. The lack of symmetry for the error estimates 
when using Virgo or Fornax as the calibrator is due to follow- 
ing: The error estimates on the slope of the regression line, 
which were propagated in the above analysis, were derived 
from the analytical expression for these terms. However, the 
errors on the slope and intercept will be under-estimated by 
the analytical expressions when the number of data points 
is small. This is seen in the FP solutions shown in Table Q 
where the jackknife and bootstrapping methods of analysis 
provide a better estimate of the true size of these errors. 
Consequently, the error estimates on the previously given 
distance moduli should be multiplied by approximately the 
ratio of the bootstrapping error (or jackknife error if this is 
larger) to the asymptotic error (Feigelson & Babu 1992). 

In addition, it should be noted that the FP slope uncer- 
tainties in Table [?] were estimated when the mixing value, 
b, was held fixed and assumed to have zero error. It should 
be remembered that it is the offset between two parallel 
planes, not two parallel lines, that is being measured. In 
order to fully estimate the uncertainties on the derived dis- 
tance modulus, one must allow for errors in the optimal FP 
slopes for all axes. In so doing, one should return to the 
bootstrapping errors obtained with the 3-dimensional data 
set (Table ^). Working in 3-dimensions, the errors associated 
with the FP exponent A, are obtained by also letting the 
exponent B (and hence 6) vary, as permitted by the boot- 
strapping analysis of the data. Therefore, the error obtained 
in the 3-dimensional analysis represents the true error esti- 
mate of the exponent A. Allowing for this, the errors to the 
solution for the Virgo-Fornax distance modulus are again 
multiplied by the ratio of the 3-D bootstrapping errors to 
the 2-D asymptotic errors. This procedure does not effect 
the optimal solution, only the associated uncertainties. The 
resulting solutions place Fornax (23±44)% (_R 1//4 ,<to), and 
(12±36)% (i? 1/n ,crtot,n) further away than Virgo. 



5 DISCUSSION 

In the interest of statistical completeness, galaxies with a 
CVD of less than 100 km s _1 were not explicitly removed 
from the analysis. Due to the greater scatter in the FP below 
this value, past studies have enforced this arbitrary cut-off in 
their galaxy sample. However, it seems plausible to attribute 
the higher level of scatter to an inadequate treatment of the 
kinetic energy for these galaxies. As this is an issue which is 
addressed in this study, it made sense to leave these galaxies 
in the analysis. 

Contributions from rotational energy to the total ki- 
netic energy of a galaxy can in principle influence the slope 
of the FP. Djorgovski & Santiago (1993) noted that gains 
from including a rotational energy measurement would be 
offset by the introduction of additional measurement errors. 
Indeed, a marginal effect has been reported many times 
(Simien & Prugniel 1992; Busarello, Longo & Feoli 1992; 
Prugniel & Simien 1994). An asymmetric trend between the 
galaxy FP residuals and their maximum rotational velocity 
has also been found by D'Onofrio et al. (1996). More quan- 
titative estimates from Prugniel & Simien (1996a,b), using 
a sample of 400 galaxies, claim that rotation can account for 
25±15% of the departure of the FP from the virial expec- 
tation. While this figure may be true for FP's constructed 
using early- type galaxies (E and SO), which includes almost 
all FP's constructed to date, this figure can be reduced by 
the exclusion of SO galaxies, known to have rotating disks. 
This was hinted at in the study by Busarello et al. (1997), 
who used only bona fide ellipticals to obtain a contribution 
of 15±28% to the tilt to the FP due to galaxy rotational 
energies. 

In this paper, the FP was constructed using a culled set 
of galaxies, removing the SOs and keeping the Es. Contribu- 
tions from possible rotational support to the kinetic energy 
term are expected to be minimal in this study of the FP. 
Four of the nine SO galaxies which were excluded did ac- 
tually have CVD's less than 100 km s _ , leaving only one 
galaxy in the remaining sample of 25 ellipticals with a CVD 
less than 100 km s _ . It may indeed be that the increased 
scatter observed for the FP at the low velocity dispersion 
end is due to the presence of SO galaxies that harbour rotat- 
ing disks with significant rotational energies that are simply 
not measured through CVD measurements. For the Virgo 
cluster, exclusion of these objects resulted in a reduction of 
some 15-35% in the scatter about the FP, as evidenced by 
the rms residuals listed in Table ti and Table a. 
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For the Virgo cluster, exclusion of the SO's caused 
the final FP solution (incorporating the observed range of 
galaxy profiles and implied dynamics) to change from i? c ,n oc 

„1.37±0. 16v-0. 76±0. 05 + „ p ^ ^1.72±0.24 v -0.74±0.09 rpy. 
a tot,n ^c.n to -tto.n OC <J totn ^ c ,„ • 1 U1S 

suggests that the inclusion of the SO population not only in- 
creases the scatter about the FP (see Tables |2| and |^) but is 
also partly responsible for the departure of the FP from the 
relation R oc <r 2 E _1 , predicted by the virial theorem. This 
is again easily understood as due to the neglect of rotational 
energies in the construction of the observed FP. It is there- 
fore suggested that in addition to treating broken structural 
and dynamical homology, the construction of the FP always 
adopt the procedure of only using bona fide elliptical galax- 
ies, unless one allows for rotational energies. 

It is found that by addressing the range of structural 
and dynamical profiles evident amongst the elliptical galaxy 
population, the departure of the FP from the plane predicted 
by the virial theorem is reduced, confirming the results in 
Paper I. In addition, the scatter of the data points about the 
FP is reduced by some 20% for the sample of Fornax ellipti- 
cals and by the same amount for the sample of Virgo E and 
SO galaxies. For the sample of Virgo ellipticals, the mean 
rms scatter about the FP is reduced by 40% after treating 
the range of galaxy structures and dynamics. This is evident 
in Tables 2, 3 and 4, which give the best fitting FP relation 
when one does and does not treat broken structural homol- 
ogy, and when one does and does not consistently measure 
the kinetic energy of the galaxies (i.e. uses a to t or CVD's, re- 
spectively). The first row of these 3 tables gives the FP which 
was constructed assuming that all galaxies are described by 
an profile and uses CVD's to represent the kinetic en- 
ergy of each galaxy; This procedure gives the greatest depar- 
tures of the FP from the virial plane. The final row in each 
Table gives the FP that was constructed allowing for a range 
of structural profiles amongst the galaxies, as described by 
the R 1/n model, and measures the kinetic energy of each 
galaxy in a consistent manner by using the computed ve- 
locity dispersion within an aperture of infinite radius; This 
procedure gives the closest agreement between the FP and 
the virial plane. In fact, the formal solution for the FP of 
the Fornax cluster (i? c , n oc cr t 2 o ? 3 n ±0 - 78 E-i 07±0 ' 30 ) is found to 
be in very close agreement with the prediction of the virial 
theorem, but does have large uncertainties associated with it 
because of the small number (7) of Fornax elliptical galaxies 
in the sample. Overall, it is concluded that addressing the 
range of structural profiles amongst the elliptical galaxies 
and better treating the galaxy dynamics reduces both the 
departure of the FP from the virial plane and the scatter 
about the FP. 

In Paper I the Sersic light profile model was fitted to a 
sample of 26 E/SO Virgo galaxies using the aperture mag- 
nitude profile data from Bower et al. (1992). Of these, 18 
galaxies are in common with the sample of 25 used here. 
Comparison between the galaxy model parameters obtained 
in each study is not straight-forward, however. The galaxy 
profile data used here is in the form of surface brightness 
measurements along the major-axis, where as Paper I uses 
circular aperture magnitudes. Consequently, the models fit- 
ted in Paper I are more heavily weighted by the inner profile, 
and by necessity incorporate the central parts of the galaxy 
light. By contrast, the treatment of the surface brightness 
profiles here gives all data points along the profile equal 



weight, and does not fit the central portion of the profile. 
In addition, the use of aperture magnitude data samples 
light from all position angles, and in so doing represents a 
mean light distribution for the galaxy as a whole. The pres- 
ence of ellipticity gradients will result in the major- and 
minor-axis light profiles being described by different shape 
parameters so that n ma j^n m i n . This is clearly evident in 
the range of shape parameters, n, for the major- and minor- 
axis, shown in the study of the Virgo ellipticals by Caon et al. 
(1993). Consequently, the mean light distribution obtained 
through the use of circular aperture magnitudes should not 
be expected to be identical to the light distribution along 
the major-axis. Finally, the data of Bower et al. (1992) was 
taken in the V-band, and that of Caon et al. (1993) in the 
B-band. Possible colour gradients, although expected to be 
small, may also contribute to differing results. 

Despite the above factors, there is still fair agreement 
between the Sersic model parameters obtained using the dif- 
ferent data sets. Galaxy scale-size and light profile shape for 
those galaxies in common between this Paper and Paper I 
are displayed in Figure |S| For comparison between the data 
sets, the Sersic model was fitted to the same radial extent 
(<40") for each data set. This outer radius cut-off comes 
from the aperture magnitude profile limit of Bower et al. 
(1992), and has been applied, for this comparison, to the 
major-axis profiles of Caon et al. (1994). The correlation 
seen between the shape parameters has a linear correlation 
coefficient of 0.55, which increases to 0.82 upon removal of 
the 3 more discrepant points that are labeled in Figure |5| 
The correlation coefficient between the effective half-light 
radii is stronger at 0.93, suggesting that the range of galaxy 
sizes, and structures, can be largely reproduced from both 
types of observational data. The galaxies which appear to 
have the more discrepant shape parameters are the flattened 
SO galaxy NGC 4550, NGC 4486, and NGC 4564 which has 
a value of n=4.64 (Paper I) lying between n ma j = 1.51 and 
n m i n — 5.08 (Caon et al. 1993) and is therefore not a dis- 
crepant point after all. 

Despite the above differences, it is interesting to com- 
pare the resulting FP's obtained using the independent data 
sets. Firstly, assuming structural homology and using R 1 ' 4 
model parameters and CVD measurements for the Virgo 
cluster of galaxies, the FP was found in Paper I to be de- 
scribed by R cA oc a i-33±o.io E -o.79±o.ii Thig ig tQ be CQm _ 
pared with R cA oc a i io±o.i4 E -o.55±o.09 found in tMg gtudy 
The exponents on the velocity dispersion term are consis- 
tent within the la joint errors, and the exponents on the 
surface brightness term are consistent within the 2er joint 
errors. Allowing for broken structural homology and con- 
sistently measuring the velocity dispersion within an aper- 
ture of infinite radius, it was found in Paper I that the FP 
for Virgo is described by J? e ,n oc o'l 44 ^ z0 ' 1 x 



should be compared with i? c , n oc ^STio.ie v -o.76±o.o 5 



^<T,4 • The 



relative consistency between the exponents is the same as 
before. It is noted that the uncertainties on the FP expo- 
nents are based upon the statistics of the linear regression 
to the data points; they do not incorporate possible mea- 
surement errors in the data points themselves. It is also 
noted that the sample of Virgo galaxies used in Paper I 
contains 7 galaxies that this study does not include, and 
this study also contains 7 galaxies not present in the sam- 
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Figure 3. The Sersic R}/ n model parameters for the effective 
radius, R e , and shape parameter, n, are displayed for 18 galax- 
ies from Paper I which are in common with with the sample in 
this study. The subscript (1) denotes that the parameters were 
obtained from Paper I, and the subscript (2) refers to the values 
obtained in this study. 



pie of Paper I. Using only those 18 galaxies in common, 
the resulting FP's were found to be in acceptable agreement 
with each other, with the FP described by the data in Pa- 

- 1.10 E -0.89±0.05 d 



per I given by the relation Rc, n oc <r ti 
that obtained using this studies data given by the relation 

o ^.1.26±0.12 V -0.81±" 
-n,c.n ^ U tnt n ^c.n 



°"toi,n • The main P° int iS that b0th 

studies show that by treating the range of luminosity and 
dynamical structures inherent amongst the early type galax- 
ies the departure of the FP from the virial plane is reduced. 

D'Onofrio et al. (1996) also constructed a FP for the 
Virgo and Fornax clusters using the data of Caon et al. 
(1994). The approach taken in this study differs in sev- 
eral ways to that of D'Onofrio et al. (1996). This work 
has used model-determined effective radii and the model- 
determined mean surface brightness within these radii, 
where as D'Onofrio et al. (1996) used model-independent 
values for these quantities (Caon et al. 1994). The model- 
determined quantities are those derived from the R 1 ' n model 



fitted to the undisturbed portion the galaxy light profile. 
This portion is free from any peculiarities of the inner light 
profiles (noted by Caon et al. 1993), is not heavily distorted 
by the presence of a disk or dust lanes, and is uncontami- 
nated from possible tidal interactions which may have per- 
turbed the original galaxy profile, particularly in the outer 
parts. The parameters used by D'Onofrio et al. (1996) are 
better estimates of the half-light radii, in the true sense of 
the term, in that they are the radii which enclose half of 
the total light from the galaxy, although the total light was 
derived by extrapolating the observational data to /is =32 
with R 1 / 4 profiles (Caon et al. 1994), The preferred dynam- 
ical term used here differs from that used by D'Onofrio et 
al. (1996). While they used CVD values, this paper has also 
used, and in fact prefers, the infinite aperture velocity dis- 
persion derived from the application of the Jeans equation 
(Ciotti 1991; Paper 1). D'Onofrio et al. (1996) gave their 
Virgo FP solution as # c cx(T ( V 26±0 ' 09 £- a70±0 - 03 . This is in 
good agreement with the FP solution derived here, which 
addressed the range of structural profiles but used CVD 
measurements. For the sample of Virgo E and SO galaxies 



the relation R c n <xa n 



was obtained here. 



D'Onofrio et al. (1996) computed a Virgo-Fornax dis- 
tance modulus of 0.50 mag (using the FP) and 0.45±0.15 
mag (using the D n -a relation). This result is in close agree- 
ment with the distance modulus derived in this study based 
upon galaxy parameters from the R 1 '* profile and using 
CVD's, where a value of 0.45±0.16 mag (asymptotic er- 
ror) was computed. Using J? 1//n model parameters and in- 
finite aperture velocity dispersion measurements, however, 
this study found the Virgo-Fornax distance modulus to be 
0.25±0.12 mag (asymptotic error). Given that D'Onofrio et 
al. (1996) allowed for different galaxy structures through us- 
ing model-independent estimates of the half-light radius and 
associated mean surface brightness, but used CVD's, the fi- 
nal estimate of 0.25 mag appears to be different because of 
the treatment of the dynamical term rather than the treat- 
ment of the structural quantities. Indeed, D'Onofrio et al. 
(1997) introduced their own aperture correction to the veloc- 
ity dispersion term, finding a reduced Virgo-Fornax distance 
modulus of 0.30±0.05 mag (formal error). 

Comparison with other methods of distance determi- 
nation yields good agreement with the distance modulus of 
+0.25±0.12 mag obtained in Section ^. The investigation 
by McMillan et al. (1993), using the planetary nebula lumi- 
nosity function, found a Virgo-Fornax distance modulus of 
+0.24T0.10 mag. Another accurate estimate has been de- 
rived using the method of surface brightness fluctuations, 
where a value of 0.20T0.08 mag has been obtained (Tonry 
et al. 1997). Using a sample of E and SO galaxies, Dressier 
et al. (1987) applied the D n -a relation to compute a value 
of +0.14±0.18 mag. An alternate approach using type la 
supernovae (Hamuy et al. 1991) found a yet lower value of 
0.09T0.14 mag, although still consistent within the errors of 
the value derived here. 

However the result presented here is in poorer agree- 
ment with studies using the globular cluster luminosity func- 
tion (-0.5±0.2 mag, Bridges, Hanes & Harris 1991) and with 
studies which have used the Tully Fisher (TF) relation for 
spiral galaxies. The H-band TF relation used by Aaronson 
et al. (1989) placed the Fornax cluster closer than Virgo, 
giving a distance modulus of -0.25±0.23 mag for Fornax rel- 
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ative to Virgo. More recently, Bureau et al. (1996) used an 
I-band TF relation. Demanding equal slopes for the TF re- 
lation in each cluster, and allowing for errors on both axes, 
they obtained a value of -0.04±0.15 mag. This result is not 
as inconsistent with the analysis in this paper as one might 
first think. It should be noted that the bootstrapping error 
analysis, performed in Section Q revealed that the estimate 
of the error based upon analytical asymptotic expressions 
was too small due to the small sample size (see Table 0). 
Correction for this increased the error estimate from 0.12 
mag to 0.27 mag. The treatment in Bureau et al. (1996), 
and indeed most papers of this kind, has similarly under- 
estimated the true size of their errors for the same reason 
- the analytical expression for the size of the errors under- 
estimates the true size of the errors when the sample size 
is small (approximately <50). In addition, it is known that 
spiral galaxies are less centrally concentrated within a clus- 
ter than the elliptical galaxies, so they may not provide as 
accurate a measurement of the cluster distance as the el- 
lipticals can. This could be particularly relevant for Virgo, 
a cluster known to have a considerable line-of-sight depth 
(Tonry, Ajhar & Luppino 1990; Jacoby et al. 1992). 

The Virgo-Fornax distance modulus obtained in this 
paper is also only marginally consistent with the study by 
Bothun, Caldwell & Schombert (1989) who used the surface 
brightness profiles of dwarf ellipticals to compute a distance 
modulus of -0.16±0.16 mag. By the nature of their study, 
they also sampled a different galaxy population than used 
here, which again may partly explain the disagreement in 
distance modulus. It is also observed that they made no 
mention of having used any numerical resampling technique 
to estimate the size of their errors. To extrapolate, this opens 
up the possibility that many past inconsistencies between 
cluster distance estimations may be simply due to an inad- 
equate, and under-estimation, of the associated errors. 

There are other effects which have been proposed to in- 
fluence the tilt of the FP. One of these is orbital anisotropy 
in the galaxies' internal velocity fields. This has been ex- 
plored in detail in the studies by Ciotti, Lanzoni & Renzini 
(1996) and Ciotti & Lanzoni (1997). Using the Jeans equa- 
tion, they performed theoretical modeling of the dynamics 
of a range of R 1 ^ 71 surface brightness profiles (n=l-10), and 
for the Hernquist (1990) and Jaffe (1983) density models. 
They explored possible variations in the orbital anisotropy 
with radius. Rejecting those models which are likely to be 
unstable, they found that variations to the tilt of the FP due 
to orbital anisotropy are confined within the known thick- 
ness of the FP and can thus have no significant influence on 
the FP tilt. 

The influence of differing stellar populations, as given 
by line strength indices such as Mg2 or broad band colours 
such as B-V, was initially determined to have little or no 
effect on the slope of the FP (Djorgovski & Davis 1987; 
Lynden-Bell et al. 1988; j0rgensen et al. 1993). However, 
more recent estimates of the influence from differing stellar 
metallicities and ages are suggesting a greater effect on the 
slope of the FP (Gregg 1992; Guzman & Lucey 1993). Prug- 
niel & Simien ( 1996a, b) combined four broad band colours 
and the Mg2 index in the extraction and computation of a 
mean measure of the influence of differing stellar popula- 
tions. They found that stellar populations may account for 
half of the observed tilt to the FP. If confirmed, this would 



easily explain any remaining tilt to the FP, and may place 
one in the unexpected position of having the FP tilt the 
other way. 



6 SUMMARY 

For a sample of Virgo and Fornax early-type galaxies, the 
surface brightness profiles from Caon et al. (1994), which 
are known to represent a range of luminosity structures, 
have been used to explore the tilt and thickness of the FP. 
This study also explores the range of luminosity structures 
present amongst the early-type galaxies using the Sersic 
Ft}/ n model rather than forcing structural homology through 
use of the R 1/4 model (see also Caon et al. 1993). The R 1/n 
light profile model has been used to parameterise the size (as 
given by the effective radius _R e ,n) and mean surface bright- 
ness (S e , n ) for each galaxy. Solving the Jeans equation, the 
associated range of dynamical structures has been computed 
and the projected velocity dispersion within an aperture of 
infinite radius (=1/3 of the virial velocity dispersion; Ciotti 
1994) has been derived. This approach measures the veloc- 
ity dispersion of each galaxy in a consistent fashion. CVD 
measurements sample different fractions of different galaxy's 
dynamical profiles, where as the infinite aperture velocity 
dispersion consistently measures galaxies in the same way, 
and therefore provides a better estimate of a galaxy's kinetic 
energy. 

The trend of increasing shape parameter n with increas- 
ing values of the model-independent half-light radius (Caon 
et al. 1993), is also shown to be present when using the R}' n 
model-determined effective radius R e ,n- A \ 2 error analysis 
of the fitted model parameters shows that this is not due 
to coupling between the model parameters implying it is a 
real physical relation. A principal component analysis re- 
veals that greater than 97% of the variance in the FP data 
set (<r,_R,E) resides within 2 dimensions, confirming the ap- 
propriateness of fitting a plane to this data. 

It is concluded that broken structural and dynamical 
homology are responsible, at least in part, for the depar- 
ture of the FP from the relation i?oc<r 2 E -1 predicted by the 
virial theorem. Having used an independent set of data ob- 
tained by different observers with a different instrument and 
reduction procedures, this confirms the results of Paper I. 
This conclusion is supported by the additional analysis of 
the Fornax cluster. 

For comparison of the the FP with the virial theorem, 
the FP was constructed using the bisector method of linear 
regression in 3 dimensions. Due to the small galaxy sam- 
ple size, a bootstrapping procedure was employed to esti- 
mate the uncertainty on the FP slope. Using a sample of 25 
Virgo E (18) and SO (7) galaxies, the FP constructed using 
i? 1 / 4 model parameters and CVD measurements was found 
to be given by the relation R cA oc .i-i°±o.i4 E -o.55±o.o9^ 
Allowing for broken structural and broken dynamical ho- 
mology, this plane was observed to change to -R e , n oc 
<7to?,n ±0 ' 16E ^n 76±0 ' 05 - The perpendicular rms scatter about 
these FP's is 0.099 to 0.077 dex respectively. 

Given that the SO galaxies likely possess significant 
rotational energy that is not dealt with through use 
of velocity dispersion measurements, they were removed 
from the sample of Virgo galaxies and the above anal- 
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ysis was rc-performed. The resulting FP was found to 
change from R eA oc ff i-i9±0-2i E -o.60±o.ii to R ^ K 

o-to[ 2 n ±0 ' 24 ^n 74±0 ' 09 . The perpendicular rms scatter to the 
plane was found to decrease by 40% from 0.084 to 0.050 with 
the fuller treatment of galaxy structure and dynamics. 

The FP for a sample of 7 elliptical galaxies in the Fornax 
cluster, was also computed under the same conditions as 
above. The Fornax FP was found to change from R aA oc 

1.69±0.83 v -0.85±0.30 . p „ 2.03±0.78 v -1.07±0.30 ' 
a Q L cA tO K , n (X <7 tot]n 2j e ,n , With 

perpendicular rms scatter of 0.064 and 0.050 respectively. 
The formal solution of the latter FP is in close agreement 
with the plane predicted by the virial theorem, and it will 
be of interest to see how/if this solution changes for a larger 
sample of galaxies. 

Constructing the Virgo FP by minimising the residuals 
of the log 7? variable and applying this to the Fornax data, 
through use of the improved Working-Hotelling confidence 
bands approach developed by Feigelson & Babu (1992), en- 
abled an estimation of the relative distance between the 
two clusters. Forcing structural homology (i.e. fitting R 1 '' 4 
profiles) and using CVD's, a Virgo-Fornax distance modu- 
lus of 0.46±0.16 mag (analytical error) was obtained. Al- 
lowing for the range of galaxy structures observed in the 
data, and using the infinite aperture velocity dispersion, 
the Virgo-Fornax distance modulus was computed to be 
0.25±0.12 mag (analytical error). A fuller estimate of the as- 
sociated uncertainties to these distance moduli was achieved 
through a numerical resampling of the 3-dimensional FP 
data set. This analysis made the Fornax cluster (23±44)% 
and (12±36)% more distant than the Virgo cluster, for the 
respective models above. 



Acknowledgments 

I thank Mauro D'Onofrio for making his galaxy profile data 
available to me in electronic form. I also thank Eric Feigel- 
son and Ron Kollgaard for providing me with their com- 
puter codes to compute the intercept offset between parallel 
data sets, and for their quick help with my queries. I am 
also grateful for the use of Eric Feigelson's computer code 
SLOPES. I thank Matthew Colless for his proof reading of 
this manuscript and his useful suggestions. This research 
has made use of the NASA/IPAC Extragalactic Database 
(NED) which is operated by the Jet Propulsion Laboratory, 
California Institute of Technology, under contract with the 
National Aeronautics and Space Administration. 



REFERENCES 

Aaronson M. ct al., 1989, ApJ, 338, 654 

Bender R., Saglia R.P., Gerhard O.E., 1994, MNRAS, 269, 785 
Bevington P.R., 1969, Data Reduction and Error Analysis for the 

Physical Sciences. New York, McGraw-Hill 
Binggeli B., Cameron L.M., 1991, A&A, 252, 27 
Bothun G.D., Caldwell N., Schombert J.M., 1989, AJ, 98, 1542 
Bower R.G., Lucey J.R., Ellis R.S., 1992, MNRAS, 254, 589 
Bridges T.J., Hanes D.A., Harris, W.E., 1991, AJ, 101, 469 
Bureau M., Mould J.R., Staveley-Smith, L., 1996, ApJ, 463, 60 
Burstein D., Heiles C, 1978, ApJ, 225, 40 
Burstein D., Heiles C, 1984, ApJS, 54, 33 

Burstein D., Davies R.L., Dressier A., Faber S.M., Stone R.P.S., 



Lynden-Bell D., Terlevich R.J., Wcgner G., 1987, ApJS, 64, 
601 

Busarello G., Longo G., Fcoli A., 1992, A&A, 262, 52 
Busarcllo G., Capaccioli M., Capozziello S., Longo G., Puddu E., 

1997, A&A, 320, 415 
Caon N., Capaccioli M., D'Onofrio M., 1993, MNRAS, 265, 1013 
Caon N., Capaccioli M., D'Onofrio M., 1994, A&AS, 106, 199 
Capaccioli M., 1989, in Corwin H G., Bottinelli L., eds, The World 

of Galaxies. Springer- Verlag, Berlin, p. 208 
Capaccioli M., Caon N., 1989, in Grosb0l P.J., Murtagh F., 

Warmels R.H. eds., Proc. ESO/ST-ECF workshop on: Data 

Analysis. ESO, Garching, p. 107 
Capaccioli M., Longo G., 1994, A&A Review, 5, 293 
Ciotti L., 1991, A&A, 249, 99 

Ciotti L., 1994, Cel. Mech. & Dynam. Astr., 60, 401 

Ciotti L., Lanzoni B., 1997, A&A, 321, 724 

Ciotti L., Lanzoni B., Rcnzini A., 1996, MNRAS, 282, 1 

Davies R.L., Efstathiou G., Fall S.M., Illingworth G., Schecter 

P.L., 1983, ApJ, 266, 41 
dc Vaucoulcurs G., 1948, Ann d'Ap., 11, 247 
de Vaucoulcurs G., 1953, MNRAS, 113, 134 

de Vaucouleurs G., dc Vaucoulcurs A., Corwin H.G.Jr., Buta R.J., 
Paturel G., Fouque P. 1991, in Third reference catalogue of 
bright galaxies. Springer- Verlag 

Djorgovski S., Davis M., 1987, ApJ, 313, 59 

Djorgovski S., Santiago B.X., 1993, in Danzigcr J. et al./, eds, 
Proc. of the ESO/EIPC Workshop on Structure, Dynamics, 
and Chemical Evolution of Early- Type Galaxies. ESO Publi- 
cation 45, p. 59 

D'Onofrio M., Capaccioli M., Caon N., 1994, MNRAS, 271, 523 
D'Onofrio M., Capaccioli M., Zaggia, S.R., Caon N., 1997, MN- 
RAS, in press 

D'Onofrio M., Longo G., Capaccioli M., 1996, in Buzzoni A., Ren- 
zini A., Serrano G., eds, ASP Conf. Ser. Vol. 86, Fresh views 
of elliptical galaxies. Astron. Soc. Pac, san Francisco, p. 143 
Dressier A., Lynden-Bell D., Burstein D., Davies R.L., Fabcr 

S.M., Terlevich R.J., Wegner G., 1987, ApJ, 313, 42 
Feigelson E.D., Babu G.J., 1992, ApJ, 397, 55 
Graham A., Colless M.M., 1997a, MNRAS, 287, 221, Paper I 
Graham A., Colless M.M., 1997b, in Arnaboldi M., Da Costa, 
G.S., Saha P., eds, ASP Conf. Ser. Vol. 116, The Second 
Stromlo Symposium - The Nature of Elliptical Galaxies. As- 
tron. Soc. Pac, san Francisco, p. 175 
Graham A., Lauer T.R., Colless M., Postman M., 1996, ApJ, 465, 
534 

Gregg M.D., 1992, ApJ, 384, 43 

Guzman R., Lucey J.R., 1993, MNRAS, 263 L47 

Hamuy M., Phillips M.M., Maza J., Wischnjewsky M., Uomoto 

A., Landolt A.U., Khatwani R., 1991, AJ, 102, 208 
Hernquist L., 1990, ApJ, 356, 359 

Illingworth G., 1983, in Athanassoula E., ed., Proc. IAU Symp. 

100, Internal Kinematics and Dynamics of Galaxies. Kluwer, 

Dordrecht, p. 257 
Isobc T., Feigelson E.D., Akritas M.G., Babu G.J., 1990, ApJ, 

364, 104 

Jacoby G.H. ct al., 1992, in ??? eds, ASP Conf. Ser. Vol. 104, ??? 

Astron. Soc. Pac, san Francisco, p. 599 
Jaffe W., 1983, MNRAS, 202, 995 

j0rgensen I., Franx M., Kjasrgaard P., 1993, ApJ, 411, 34 
Lynden-Bell D., Burstein D., Davies R.L., Dressier A., Fabcr 
S.M., Terlevich R., Wegner G., 1988a, in Pritchet C.J., Van- 
derBergh S., eds, ASP Conf. Proc. 4, The Extragalactic Dis- 
tance Scale, ASP, San Francisco, p. 307 
McElroy D.B., 1995, ApJS, 100, 105 

McMillan R., Ciardullo R., Jacoby G.H., 1993, ApJ, 416, 62 
Michard R., 1985, A&AS, 59, 205 

Murtagh F., Heck A., 1987, Multivariate Data Analysis. Reidel, 
Dordrecht 



© 0000 RAS, MNRAS 000, 000-000 



A Virgo-Fornax distance modulus 

Prugniel Ph., Simien F., 1994, A&A, 281, LI 

Pmgniel Ph., Simien F., 1996a, A&A, 309, 749 

Prugniel Ph., Simien F., 1996b, in Buzzoni A., Rcnzini A., Ser- 
rano G., cds, ASP Conf. Ser. Vol. 86, Fresh views of elliptical 
galaxies. Astron. Soc. Pac., san Francisco, p. 143 

Schombert J.M., 1986, ApJS, 60, 603 

Sersic J.-L., 1968, Atlas de Galaxias Australes, Observatorio As- 

tronomico, Cordoba 
Simien F., Prugniel Ph., 1992, in Longo G., Capaccioli M., 

Busarcllo G., eds, Morphological and physical classification 

of galaxies. Kluwcr, Dordrecht, p. 431 
Tonry J.L., Ajhar E.A., Luppino G.A., 1990, AJ, 100, 1416 
Tonry J.L., Blakeslee J. P., Ajhar E.A., Dressier A., 1997, ApJ, 

475, 399 

Working H., Hotelling H., 1929, J. Am. Stat. Assoc. Suppl., 24 73 



© 0000 RAS, MNRAS 000, 000-000 



